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It is the purpose of the present article to show that so-called network models, originally designed 
to describe static properties of disordered electronic systems, can be easily generalized to quantum- 
dynamical models, which then allow for an investigation of dynamical and spectral aspects. This 
concept is exemplified by the Chalker-Coddington model for the Quantum Hall effect and a three- 
dimensional generalization of it. We simulate phase coherent diffusion of wave packets and consider 
spatial and spectral correlations of network eigenstates as well as the distribution of (quasi-) energy 
levels. Apart from that it is demonstrated how network models can be used to determine two-point 
conductances. Our numerical calculations for the three-dimensional model at the Metal-Insulator 
transition point delivers among others an anomalous diffusion exponent of 77 = 3 — D2 = 1.7 ± 0.1. 
The methods presented here in detail have been used partially in earlier work. 
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I. INTRODUCTION 

Models formulated in scattering theoretical terms have 
been proven to be convenient for the investigation of 
quantum particle- or classical wave-propagation in dis- 
ordered media. Examples related to Anderson localiza- 
tion PI are the scattering model of a one-dimensional 
(ID) disordered conductor used by Anderson ct al. M, 
it's three-dimensional (3D) version introduced by Shapiro 
PI , and the Chalker-Coddington network of the Quantum 
Hall transition B. Numerous generalizations and varia- 
tions of these models have been used in order to study 
various topics in different fields p|-p^. 

The common concept in these approaches is the scat- 
tering theoretical formulation of the problem. This leads 
to models which can be described as networks consisting 
of an array of local scattering centers, linked together by 
ID channels [g|. Usually, the static properties of these 
models are investigated, like probabilities for scattering 
through the network or correlations of stationary scat- 
tering states. 

It is the purpose of the present article to show that in 
general these models can be easily extended to dynami- 
cal network models, which then allow to study dynami- 
cal as well as spectral aspects of the respective problem 
under consideration. In particular, we consider dynam- 
ical versions of the Chalker-Coddington model and of a 
3D generalization of it Iin] . Within these models we in- 
vestigate numerically phase coherent diffusion of wave- 
packets, whereby special attention is paid to anomalous 
diffusion near the localization-delocalization transition 
point. It is demonstrated how spectral information can 
be extracted from these models, like (quasi-)energy spec- 
tra or the local density of states. Moreover, we present 
an example within the framework of a network model 
for determining conductances of an integer quantum Hall 
system contacted by point-contacts. 



A network model is mathematically formulated by an 
explicitly defined unitary matrix, the network operator 
U, that determines the scattering processes between in- 
ternal network channels. By definition, in the dynam- 
ical model this operator U will describe the evolution 
of network-states ^ over a microscopic time interval r, 
'i' — > U'i'. By this, a discrete time evolution is given: 
*(fcr) = U'"^{0), where k ^ 0,1,2,.... Further, net- 
work eigenstates $„ are defined as eigenvectors of U. 
The phases aj„ G [0, 27r[ of the unimodular eigenvalues 
exp(jix'„) will be interpreted as quasi-energy levels. 

The main part of this work consists of numerical sim- 
ulations within the two aforementioned models in order 
to demonstrate that the concept of dynamical networks 
is self-consistent and, above all, that it delivers results 
consistent with previous results and general theory of 
disordered systems. These simulations also clearly show 
that network models are convenient for numerical pur- 
poses. They allow e.g. to calculate efficiently the (fully 
phase coherent) time evolution of wave-packets diffusing 
in comparably large systems. Furthermore, the parti- 
cle energy (to be distinguished from the network specific 
quasi-energy), becomes simply a parameter of the model. 
Due to this it is possible to focus at specific and sharply 
defined energy regimes, which will turn out to be advan- 
tageous when investigating critical level statistics. 

Our approach is related to that of Edrei et al. [g| , where 
a network model has been used for calculating wave prop- 
agation through random media. Actually, the definition 
of network states and operator used here are, in princi- 
ple, identical to those in |q|. However, their work concen- 
trated on stationary states of networks with open bound- 
aries in order to determine transmission coefficients. In 
contrast to that, here we study the time evolution explic- 
itly and consider in particular closed systems. The latter 
offers the opportunity to define and to investigate corre- 
lations of network eigenstates and quasi-energy levels. 



Viewed as a time evolution operator, the network oper- 
ator is also related to so called quantum maps or Floquet 
operators, which attracted recently considerable atten- 
tion in quantum chaology ||ll|. It seems to be that dy- 
namical networks and quantum maps are identical con- 
cepts, applied in different physical contexts: the former 
to diffusive, the latter to chaotic systems. However, the 
precise relation between these two classes of systems is 
not entirely clear to us yet and deserves further investi- 
gation in the future. 

The concepts presented here have been applied par- 
tially in previous work, e.g. for determining critical eigen- 
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, local density of states |l^, c ritical level statis- 



|15| and quantum diffusion |16| . Here we explain 
in more detail the numerical methods used therein. How- 
ever, also the reader not familiar with network models 
might find this paper worth reading, since it deals with 
current topics in the field of critical disordered systems, 
like the multifractality of the local density of states, crit- 
ical diffusion or energy level correlations. 

The article is organized as follows: We start with a def- 
inition of dynamical network models, which tries to be a 
compromise between exactness and generality on the one 
hand and the amount of formal effort on the other. Next 
is the introduction of the two specific models mentioned 
above, within which several numerical simulations will 
be performed. Beginning with the simulation of diffusing 
wave packets we come to the calculation of the local spec- 
tral density. After that we turn to the distribution of net- 
work quasi-energy levels. The next section is concerned 
with correlations in the local density of critical network 
states, followed by a section in which the network model 
is used for the calculation of two-point conductances. We 
conclude with a summary and general remarks. 



II. GENERAL DEFINITIONS 

Consider a d-dimensional array of scattering centers, 
where each of them has n outgoing and incoming chan- 
nels. Let them be labeled by a lattice index j and char- 
acterized hy nxn unitary matrices Sj. Then, connecting 
in a certain manner outgoing with incoming channels of 
neighboring scatterers, one ends up with a network (Fig. 
H) . We denote its scattering centers as nodes and internal 
channels as links, labeled by another index I. 

A network is closed when all channels start from and 
end in network nodes, otherwise it is open. In the latter 
case some channels start from or end somewhere outside 
the network, forming input (/,«) and output channels 
{Oi) of the network. Thereby, the number Nc of input 
channels always equals the number of output channels. 
Such an open network can be viewed as a scatterer with a 
complex internal structure and the transmission through 
it can be described by an Nc dimensional matrix S, which 
maps the input amplitudes in /„ to the output channels 
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FIG. 1. The basic module of a network are elementary scat- 
terers of the same type. 



Oi . Both types of networks are illustrated by examples 
shown in Fig. g. 




FIG. 2. Several networks of various dimensions and channel 
numbers n: (a) a ID network with n — 2 modeling a ID chain 
of scatterers. (b) the two dimensional Chalker-Coddington 
network for the integer quantum Hall Effect, n = 2. (c) the 
same as (b) , but with spin degree of freedom, thus here n = 4. 
The network (a) is open, (b) shows an example of a closed 
network with reflective boundaries whereas the closed network 
(c) has periodic boundary conditions. 

The open network (a) would be for example appro- 
priate for the study of transmission or conductance [|6|, 
whereas the closed networks (b) and (c) could be used 
e.g. for determining eigenstates, as we will see in the 
following sections. 

One can choose refiecting or periodic boundaries for 
the network. Reflecting boundaries are realized by bend- 
ing the outgoing channels of one side to neighbored in- 
coming ones of the same side, periodic boundaries by 
linking the outgoing channels to the incoming ones of 
the opposite side (see Fig. g, (b) and (c) ). 

A network state can be introduced as the set of com- 
plex amplitudes ipi for all links I — 1, . . . ,N, '^ — {'ipi}. 
Per definition, the network states form a vector space, 
spanned e.g. by the orthonormal unit states e;/ = {^I'l}- 
The scalar product is given by (<!', ^) = J^'t'l^i- When 
the states are interpreted as wave functions of a particle, 
the squared amplitude \ipi\'^ — |(e;,4')p of a normalized 
state denotes the probability to find the particle on link 
I. 

In physically realized networks only those states \l/ are 



feasible, which have Unk ampHtudes tpi obeying the scat- 
tering condition at every node of the network. This 
means that for every outgoing hnk I of a node j with 
a scattering matrix Sj = {tj„^} the equation 
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(1) 



has to be fulfilled, whereby m runs over the n incoming 
links of that scatterer j from which link I exits (indicated 
by the symbol m : I). Stationary or steady states of the 
network are states satisfying this condition at each node. 
The right hand side of the stationarity condition Eq. 
(|l|) can be read as the action of an operator U on the 
state * = {^,}, 

N 
k—1 m:l 

Notice that the second sum contains only n summands, 
and not N (number of all links in the network) as in the 
first. This implies that only n of of the N coefficients 
{Uik}k=i,....N are non-vanishing (i.e. Uik is a sparse N- 
dimensional matrix). For closed networks, this operator 
U is unitary due to the unitarity of the local scattering 
matrices Sj. Obviously, the stationarity condition for a 
state ^ can then be written as 



t/* = *. 



(2) 



For open networks U is no longer unitary, because in 
that case the probability density conservation is violated 
by the leakage at the outgoing channels. In this case U 
can still be used to formulate the stationarity condition, 
however, additional care has to be taken of the open input 
and output channels. We will discuss this in detail for 



one example in Sec. VI] 



It is instructive to consider the action of [/ on a basic 
state e/, 
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where the sum runs over all outgoing links m of scatterer 
j in which I enters (denoted by / : m) . Operator U simply 
maps the incoming amplitude to the outgoing links of 
the scatterer according to the transmission coefficients 
^Im (Fig- 0)- A corresponding physical process can be 
thought of as the scattering of an incident wave-packet 
at a node into the outgoing channels after a characteristic 
time r. Therefore it is suggestive to consider the action 
of U on an arbitrary network state ^ as its time evolution 
over a certain microscopic time interval r. 



*(to + T) = C/*(to). 



(4) 



By this, a dynamic for the network is defined by a step- 
wise time evolution 




FIG. 3. The action of [/ on a basic state e;. It maps the 
incoming amphtude into the outgoing channels m, n accord- 
ing to the transmission coefBcients tmi , t„i . The correspond- 
ing physical process can be thought of as the scattering of 
an incident wave packet into the outgoing channels after a 
characteristic time r. 



*(to + fcr) - C/'=*(io) 



(5) 



for an integer number of time steps k = 0, ±1, ±2, . . .. 
(In the following we often set r = 1 and then write 
\l/(i) = f/*\E'(0).) Apparently, stationary states remain 
unchanged by these dynamics, as it should be. 

According to these dynamics, for closed networks a real 
spectrum {aj„} is given by the phases w„ of the complex 
eigenvalues of U, 



[/$„ = e*""$„. 



(6) 



In order to avoid ambiguity we restrict the phases or 
quasi-energies to [0, 2tt[. With respect to the dynamics 
(^ the corresponding eigenvectors $„ oscillate with fre- 
quencies UJn, 



$„(i) = [/*$„(0) - e'"*$„(0). 



(7) 



Additionally, we define a local density of network states 
(LDOS) by 

Pi(u;) = ^ |0„_;pJ(w - w„). (8) 



III. SPECIAL NETWORK MODELS 

After we have introduced network models in general, 
let us now consider two special models in more detail: 
the two-dimensional (2D) network introduced by Chalker 
and Coddington Q to describe Quantum-Hall systems 
and a 3D network very similar to the model proposed 
by Chalker and Dohmen to describe multi-layer systems 
|lO| . Originally, both models are purely static, but they 
can be easily extended to dynamical models by providing 
them with the discrete time evolution (m, as we will do 
here. 

The Chalker-Coddington network is designed in the 
style of the networks introduced by Anderson et al. H 
and Shapiro 0] in the early 1980s, where a scattering 
formalism is used to describe Anderson localization. It 
provides a semi-classical description of a 2D electron in 
a strong perpendicular magnetic field B and a potential 



V. The random potential V is assumed to be smooth 
with a correlation length a large compared to the mag- 
netic length Ic |17|-p^ . Classically, the electron executes 
a fast cyclotron motion on a circle of radius Zc, the cen- 
ter of which drifts slowly along contours V{r) — const.. 
For Ic "^ a it is justified to separate cyclotron and drift 
motion. Quantizing the first yields the Landau energies 
En = hu;{n + ^). The drift of the center coordinate 
can be treated semi-classically by a WKB- or stationary- 
phase-approximation ||l^.ll^,pO] , leading to the following 
picture: The probability density of an eigenstate at en- 
ergy E in the nth Landau band is concentrated along 
contours ^(q) = E — En of the disorder potential. In 
local coordinates u, v parallel (m) and transverse (v) to a 
contour I the wave function 'I'£;(r) can be well approxi- 
mated by 



C/(S„)*„ = ^r. 



(U) 



^_B(u,-y) = Vi 



\/vd{u) 



Mv/Qe'^^^'^\ 



(9) 



where ipnix) denotes the nth Hermite polynomial and 
(j){u, v) a gauge dependent phase. Because of the nor- 
malization by the local drift velocity Vd{u), the squared 
modulus of the complex coefficient V'/ gives the net cur- 
rent transported by the state along the contour I in di- 
rection of the drift motion. The approximation is valid 
everywhere where the contour line is well separated from 
others by a distance large compared to Ic. At saddle 
points of the disorder potential lying at energies close 
to £^ — En, such that the distance between two con- 
tour lines becomes of order l^ it breaks down. However, 
at these points the tunneling between different contour 
lines can be described by local 2x2 scattering matrices 
Sj{E) ~ {tl^i{E)}, relating the coefficients of incoming 
and outgoing currents by 
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(10) 



at all saddle points. For a given potential the scatter- 
ing coefficients for the saddle points can be determined 
as functions of the energy E by semi-classical methods, 
as has been shown and explicitly worked out by Fertig 
pO| . Further, knowing the matrices S{E), the scattering 
conditions ( [lO| ) form a closed system of equations for the 
coefficients ipi from which eigenenergies En and eigen- 
state coefficients -0; can be determined |20[| . 

Now, identifying the saddle points with nodes and the 
approximated contour states (0) with links, one arrives at 
a properly defined network. Finite isolated QH-Systems 
correspond to finite closed networks, systems coupled to 
external leads are described by open networks. Network 
states ^ = {ipi} determine wave functions of the real 
system via the expressions M). Moreover, the scattering 
condition (nfl) fulfilled by eigenstates of the QH-system 
is just the stationarity condition of the network; thus 
eigenenergies En and -states ^„ are determined by 



corresponding to Eq. (4.7) in Ref. ^^. Here the network 
operator U{E) is energy dependent, since the coefficients 
are functions of energy. 

The Chalker-Coddington model is basically identical 
to such a network, but exhibits some simplifications: It's 
nodes and links form a regular square lattice, in contrast 
to saddle points and contour lines of a real smooth disor- 
der potential, which may give a more irregular network. 
Further, for all nodes the transmission amplitudes for 
the scattering into the left outgoing links is set to a con- 
stant value r+, and for the scattering into the right link 
to r_ = 1 — T+. The scattering phases ipim = argt;™ 
are random variables, homogeneously distributed over 
[0, 27r], as in the case of the network model for a real 
disorder potential. The scattering amplitude T_|_ acts as 
model parameter. 

With these simplifications the network can be viewed 
as one corresponding to a QH-system with an egg-carton- 
like disorder potential. The basic length a of this po- 
tential is assumed to be large compared to Ic. Also it 
should exhibit small deviations from regularity, such that 
the flux penetrating the plaquettes fluctuates strongly on 
scales of $0j leading to random scattering phases cj)im. 
According to this picture we parameterize the transmis- 
sion amplitude by T± ^ {1 + exp{±E/Et), where E de- 
notes the electron energy (minus the cyclotron energy 
hwc{n + 1/2)) and Et a characteristic tunneling energy 
of the saddle points. Thus, the Chalker-Coddington net- 
work is build up by nodes with scattering matrices 



s,= 











t+{e) t_(e) 
-t_(e) t+{e) 



(12) 



where t±{e) = (1 -t-e"*"^)^^/^ and s = E/Et. These nodes 
are arranged on a regular square lattice with orientation 
as shown in Fig. 0(b). For fixed e the scattering phases 
ipj{e),Lp'Ae) are considered to be random numbers ho- 
mogeneously distributed over [0,27r]; except for this they 
are not further specified. As it has been shown by the 
originators, the network exhibits delocalized states only 
at a singular critical point determined by |t_| = |i+|, i.e. 
at Ec = 0. For e ^ the eigenstates have a finite lo- 
calization length S,, diverging with a universal exponent 
V « 2.3 when e approaches Ec = 0, ^ oc |e|~'' [Q. 

The e-dependent network operator of the model is 
given by the coefficients of the scattering matrix (|lj) via 
Eq. (H). Then the network dynamics are 



^(i) = C/*(e)^(0), i==0,±l,±2. 



(13) 



e-dependent quasi-energies ct;„(e) and eigenstates $«(£) 
are defined by 



f/(e)<i>„(e)=e^""(^)<i>„(e), 
the local density of network states is 



(14) 



pS.uj) 






|0„,i(e)p(5(w-w„(£)). 



(15) 



Quenched disorder is represented by the random 
phases 0j. Thus the disorder average of a quantity A 
is given by the integral over aU random-phase configu- 
rations, [A\ = /o'^ Hi "^^i*^™}- ^^ usual, in numer- 
ical calculations this ensemble average is approximated 
by the average over a finite set of disorder configura- 
tions {(pm}- In case of local observables Ai the statis- 
tics are often improved by taking also a system average 

A generalized version of this model is obtained, when 
also the transmission amplitudes become random vari- 
ables, e.g. by t-'^ = {1 + exp(e — fj,j))~^/'^, where ^j G 
[—W/2, W/2] is random and independent for each node 
j. This type of disorder corresponds to fluctuating saddle 
point energies of the random potential, which brings in 
aspects of classical percolation |7|,p|j26| . 

The 3D model we will investigate is very similar to that 
of Chalker and Dohmen ||l^ . It consists of a stack of the 
2D Chalker-Coddington networks we just have consid- 
ered, now coupled into the vertical direction by interlayer 
links as shown in Fig. Q. However, the local scattering 




FIG. 4. Structure of the 3D network model. 

matrices are defined as in the 2D network; so here the 
scattering amplitudes are also determined by the energy- 
parameter £ and the disorder is represented by random 
scattering phases. Due to the inter-layer coupling, the 
singular delocalized point aX E = oi the 2D network 
becomes an extended metallic band ranging from —Ec to 
Er, Er > 0. Outside this band the states are localized 



eigenvalues En and eigenvectors (/>„ of the Hamiltonian. 
After that, the state at any time i of a particle initially 
in state '0(0) can be calculated simply by 



V'(i)=5]e-^"*(0„,^(O)). 



(16) 



In this procedure, the diagonalization of the Hamiltonian 
takes the main part of the numerical effort and limits the 
application of this method to rather small system sizes [|. 
The dynamical network model offers a far better way to 
simulate quantum mechanical time evolution. Using the 
network dynamics Eq.(P), one can evolve any arbitrary 
state \l/(0) in discrete time steps by iterative multiplica- 
tion of ^(0) with the network operator U, 



*(i) = t/**(0) = L/*(i-l), t^ 0,1,2, 



(17) 



Since multiplication of vectors by matrices is one of the 
simplest tasks computers can perform (particularly when 
most of the matrix entries vanish, as in our case) the 
dynamics are very suitable for numerical applications. In 
fact, by this the simulation of the quantum mechanical 
time evolution becomes as easy as the simulation of a 
classical system with discrete dynamics. 

Let us consider numerical calculations for the 2D and 
3D networks described in the last section. In both cases 
we investigate the spreading of network states initially 
confined to a few links in the center of a squared (cu- 
bic) network of volume L''. The diffusion of such states 
is described by the disorder averaged ([...]) probabihty 
density 



p{r,t) ^ [\^ir,t)\' 



(18) 



{ip(r,t) understood as the link amplitude ipi(t) at link I 
closest to coordinate r.) We will determine this density 
for states at energies in the metallic, the insulating and 
in the transition regime. The latter regime is especially 
interesting, since at this energy interference effects start 
to destroy the diffusion m. 

Before the setups are explained in detail and results 
are presented, let us briefly summarize some results con- 
cerning quantum diffusion in disordered systems. For 
the transition point, these results are gained mainly from 
scaling theory ||8|,|9). 

At energies deep in the metallic regime the density 
p{r, t) is given by the Gaussian density. 



IV. PHASE COHERENT DIFFUSION IN 
DISORDERED SYSTEMS 

The standard method to determine numerically a 
quantum particle's time evolution is first to calculate all 



* For example, the Hamiltonian for a 3D system oi LxLxL = 
50 X 50 X 50 lattice points possesses 50^ energy levels and 
as many eigenvectors with 50 complex entries each. This 
yields approximately (50^)^ ~ 1.6 x 10^'' complex numbers, 
equivalent to at least about 100 Giga Byte, which to calculate 
or to back up is a hard job, also for todays computers. 



p(r, t) 



{4TrDty/2 



^-r^/iDt 



(19) 



where D is the diffusion constant (ordinary diffusion) 
Accordingly, the gth moment of the density grows like 



Mq{t) = / d'^rr'^p{r,t)oit'"^ 



(20) 



with exponent h = 1/2. 

In the vicinity of the transition point, scaUng theory 
predicts an algebraic decay of the density at short dis- 
tances r <C {t/pcY/'^ 



p{r,t)oct-°'/'^r°'-'^, 



(21) 



where D2 is the correlation dimension of critical eigen- 
states and pc the density of states at the transition point 
[|^-|9|. At large distances, r > {tjp^)'^!'^, the decrease 
of the density is still exponential. ^From this follows that 
at the critical point the gth moment scales like [p9[ 



Af,(t)oci«/'*. 



(22) 



Therefore, in d > 2 dimensions, the exponent h shrinks 
down from 1/2 in the metallic regime to the lower crit- 
ical value l/d when the transition point is approached, 
indicating the transition to the localized regime, where 
diffusion ceases. 

In the localized regime, the wave-packet spreads out 
until its size has become of the order of the localization 
length ^. Up to this point the density p(r, t) evolves in the 
same way as at the critical point. At larger times, t ^ 
Pc^^ I the density equilibrates in a volume ^'', centered 
around the starting point. 

Now, let us consider the outcomes of our computer 
simulations. We begin with the diffusion on the Chalker- 
Coddington network at critical energy E^ = 0. The ini- 
tial state VE'(O) is located on four links labeled by 0,1,2,3 
in the center of a squared closed network with periodic 
boundary conditions, 



1 



*(0) ^ -(eo + 61+62+63) 



(23) 



This state is subjected to the time evolution Eq. (P), 
whereby E — Ec- Figure visualizes the diffusion of 
a particle on a 40 x 40 node network and also demon- 
strates the numerical precision of the simulations. The 



FIG. 5. A state evolving backwards after it has moved 350 
time steps forward. The initial state is reconstructed to a 
high accuracy (last picture of the lower row). 

first picture in the upper row shows the density [^^(r, t)p 
after 350 time steps. The initial state ^(0) has been 
almost equilibrated over the whole lattice. In order to 



test that this equilibration is by no means influenced or 
even caused by accumulating numerical errors, we ap- 
plied successively the adjoint operator [/+ — U~^ on the 
final state ^(350). Doing this, the state evolves back (see 
Fig. 0) and returns after another 350 steps to the initial 
state with high accuracy, |*' - *| < lO"^'^. [] 

Fig. shows the gth moment Mq{t) to the power of 
2/q for q =2,4 and 6 at critical energy Ec = and in the 
strong localized regime at energy E = 3Et. The data for 
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FIG. 6. The 2nd (O), 4th (+) and 6th (D) moment (to 
the power of 2/2 2/4 and 2/6, respectively) of spreading 
wave-packets at the transition point at energy E = Ec — 
and in the strong localized regime, E — 3Et (where A, *, x 
correspond to the 2nd, 4th and 6th moment). The dashed 
lines have slope 1. 

Ec are gained from 5 networks of 300x300 nodes with 
periodic boundary condition. Since the moments are ob- 
tained by integration over the entire system, their fluc- 
tuations with respect to different disorder configurations 
are quite small; the average of 5 configurations yields mo- 
ments proportional to Mq{t) ex i'^'', with an exponent h 
very close to l/d = 1/2 (see Fig. ^). At times greater 
than the diffusion time of order Ini^ = ln(L/2)^ « 9 
the increase goes down, since then the extension of the 
wave-packet has reached the system size L. Away from 
the critical point, the moments saturate quickly at cer- 



^The stability of the quantum mechanical time evolution 
against numerical errors seems to be quite astonishing, since it 
is in sharp contrast to the numerical instability of the classical 
counterpart. Because of the chaotic dynamics, in the latter 
the numerical uncertainty in position and impulse makes the 
trajectory of a classical particle completely indefinite after a 
few time steps. However, the dynamical stability of quantum- 
chaotic systems is well-known and has been first observed by 
Shepelyansky [30]. 



tain values, reflecting the absence of diffusion at length 
scales greater than the localization length. 

Next, we consider the probability of the particle to 
return to the initial volume formed by the links I ~ 
0,1,2,3, 



Pit) = E li'i 



(24) 



;=o 



Since this return probability fluctuates strongly with re- 
spect to disorder, it is necessary to average over many 
configurations. The curve plotted in Fig. results from 
400 different critical networks of 100 x 100 nodes each. 
The return probability decays with a power law, p{t) oc 
i^, with exponent x = 0.76 ± 0.01. 




FIG. 7. The averaged return probability p{t) taken from 
400 networks of 100 x 100 nodes each. The data fits to a 
straight line with slope x — 0.76 ± 0.01. The dashed straight 
line with slope -1 corresponds to normal diflusion. 



In accordance to Eq. (|21j), this exponent fits well to 
the correlation dimension for critical states, which is ex- 
actly twice that, D2/2 = 0.75 ± 0.025 (see Sec. Wf). 



The straight line with slope -1 corresponds to the decay 
of the return probability one would expect for classical 
diffusion, governed by Eq. (|l9|). 

Finally, we discuss the density p(r, t) for fixed times 
t as a function of r. In order to reduce the numerical 
effort, we average p{r, t) over a ring of radius r around 
the origin: 

Pt{r)^— p{\r\,ip,t)dip. (25) 

Taking for granted that the disorder averaged density 
p(\r\,Lp,t) is radially symmetric, the new density pt{r) 
exhibits the same scaling behavior as p{r,t). However, 
due to the integration over the ring pt{r) fluctuates much 
less than p{r, t) , which considerably reduces the number 



of configurations necessary for a satisfyingly small sta- 
tistical error. Because of Eq. (pi[), the curves pt{r) for 
different times t should coincide for small radii r after 
rescaling them by 



Ptir) 



D2/2 



Ptir). 



(26) 



The averaged density ptir) of five critical systems of 
300x300 nodes and times t ranging from t = 1500 up 
to 15000 steps shows that this is indeed the case (Fig. ||). 
The algebraic decay of the density for small r changes 




FIG. 8. The rescaled density t"^'^p{r,t) for times t = 1500 
(O) up to 15000 (*). The straight lines corresponds to 



to an exponential decay at radii larger than a crossover 
length, which increases with time t. The data plotted in 
Fig. for small r fits to r^^"'', with d — D2 ^ 0.5, as 
predicted by Eq. (^. 

For the 3D model we carried out analogous simula- 
tions. We first determined the moments Mq{t) for states 
initially restricted to a plaquette in the center of a cu- 
bic network of L^ — 50'^ nodes. In Fig. ^ the mo- 
ments Mq{t)'^/'^ are plotted for q = 2,4, 6 at energies 
E — 1.5,3.4 and 6, measured in units of the tunnel en- 
ergy Et- The average of 5 disorder configurations has 
been taken. At energy E — l.b the moments scale with 
an exponent h = 1/2, indicating pure metallic diffusion 
according to Eq. (n9|) , for E = Q the moments converge to 
finite values small compared to L', which shows clearly 
the absence of diffusion, i.e. localization. For energies 
E = 3.4 ± 0.1 we found an anomalous reduced expo- 
nent h very close to 1/3. Consequently, energy E — 1.5 
must be in the metallic and .E = 6 in the strongly lo- 
calized regime. The transition energy we determined to 
be Ec — 3.4 ± 0.1, since at energies 3.3 and 3.5 signifi- 
cant deviations from the critical scaling Mq{t) ex f^"^ are 
already visible. 

The return probability p{t) oc r~^ shown in Fig. 10 re- 
sults from 200 systems of size L^ = 50'^ at critical energy 
E = 3.4. From these data we obtained x = 0.43 ± 0.04 




in the localized regime shows a strong, exponential de- 
cay. Both densities do not evolve further in time. At 
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Int 
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FIG. 9. The second moment M2{t) at energies E = 1.5 
(filled □), 3.4 (□) and 6 (A), corresponding to the metal- 
lic, transition and isolating regime of the 3D network model, 
respectively. 




FIG. 10. The critical return probability decays propor- 
tional ex t-^s/^, D2/3 = 0.43±0.04. Data obtained from 200 
systems of 50^ nodes each. The dotted line with slope —3/2 
corresponds to normal diffusion. 



and thus D2 — 3x — 1.3 ± 0.12 according to Eq. (21). 
For the 3d Anderson transition this correlation dimen- 
sion has not yet been determined very precisely. The val- 
ues obtained in previous publications are D2 = 1.7 ± 0.3 
|l|, £"2 ~ 1.45 - 1.8 m and D2 = 1.7 ±0.2 ||. An 
analytical calculation using an epsilon expansion yields 
D2 = 2-6 = 2- (3-2) = 1 m. It should be noted that 
the network model corresponds to a system with a mag- 
netic field while in Ref. JS^JS^,^ time reversal systems 
in the absence of magnetic fields were considered. 

Finally, Fig. |l^ shows the density profile p{r, t) at time 
t = 20000 of states initially located in the center of the 
network in the metallic, critical and localized regime. At 
this time, metallic states at E = I have been completely 
equilibrated over the entire system, seen by the constant 
density p{r,t), whereas the density of states at E = 6 
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FIG. 11. The density p(r,t) at time t = 20000 of states 
initially located in the center of a 3D network of 50 nodes. 
The curves correspond to energies E = lEt (metallic, □), 
Ec — 3AEt (critical, filled □) and E = 6Et (localized regime, 

A). 

the transition point, just at time t = 20000 the bulk 
of the spreading wave-packets cover the whole lattice. 
Therefore, the corresponding density decays algebraically 
oc r~^ according to Eq. ( pl| ) with exponent y close to 
rj^3-D2 = 1.7 ±0.1. 

To summarize the outcome of the simulations, both the 
2D and 3D dynamical network models show precisely the 
expected diffusive behavior. The metallic regime of the 
3D model shows normal diffusion with a Gaussian density 
p{r,t). At the critical point in both systems the simula- 
tions fit very well to the anomalous diffusion predicted by 
scaling theory. It is also quite remarkable that the sim- 
ple, locally defined dynamics (H) lead immediately to the 
well-known but highly non-trivial Anderson-localization, 
caused by long ranged destructive interferences. 

Since dynamical network models correctly describe the 
dynamics of disordered systems in the metallic, local- 
ized and even in the transition regime, their generat- 
ing operator U must necessarily contain relevant infor- 
mation on the spectral properties in the corresponding 
regimes. Certainly, this information is encoded in the 
quasi-spectrum {uj„} of U{E), on which we will focus 
now. 



V. SPECTRAL STATISTICS 

We start with the investigation of the local spectral 
density Po{lo) of the two models under consideration. Af- 
ter that we turn over to the statistics of quasi-energy- level 
distribution in the Chalker-Coddington network. As be- 



fore, particular attention is paid to the metal-insulator 
transition point. 

A. Local Spectral Density 

The local spectral density can be calculated efficiently 
by iteration in nearly the same way as before the re- 
turn probability p(t). Instead of the return probability 
Eq. (|2J) one has to determine the complex amplitude 
ip{0,t) = (eo,C^*eo), i.e. the propagator Goo{t)- In prin- 
ciple, Fourier-transforming the latter immediately yields 
the spectral density po{^) = J2n \^n{0)\'^S{iLj — oJn) be- 
cause of 

Goo(t) = (eo, C/*eo) = ^ e'""*|$„(0)p (27) 

duj e"^'po{cj). (28) 



Thus, to obtain the local spectral density po{i-o) we gen- 
erate the sequence Goo{t),t = 0, 1, ... T successively by 
applying U on the state Gq. Before doing the numeri- 
cal Fourier transformation, we softened the sharp edges 
of the time-window from to T by multiplication with 
1/2 — l/2cos(27rt/r) (the so called Hamming window), 
in order to reduce artifacts due to the finite time range. 
Additionally, to ensure numerical stability, we damped 
the Green function with the exponential e~*''^, whereby 
the damping time r fulfills 1 <C t ^ T. Then, perform- 
ing a discrete Fourier transformation on this windowed 
and damped Green function and taking the real part of 
the result we get the quantity 



P'q{uj) = y^ |(/)„(0)p(5^(a; - UJn), UJn = 



2-K 



n, 



(29) 



where Sr{uj) is a peaked function with maximum value 
at (jj = and width of order l/r. For sufficiently large 
times r and T ^ t, this quantity converges to the wanted 
density po{uj). 

By this method, the spectral densities plotted in 
Fig. |2|of a Chalker-Coddington network of 50 x 50 nodes 
have been calculated, whereby we chose r = 2000 and 
T = 32768. In the upper spectrum, at energy E = 2Et, 
the density is almost concentrated at only a few ener- 
gies, reflecting strong localization. The lower, critical 
spectrum at ii^c = fluctuates strongly and turns out to 
be scale invariant, as we will see in a moment. 

The same method applied to the 3D network model 
considered above gave the three spectra shown in Fig. \[q , 
corresponding to the localized, critical and delocalized 
regime. For the calculations we used systems of 50 x 50 x 
50 nodes, a maximum time T = 32768 and t = 4000. At 
energy £■ = 6 in the localized regime and at the transi- 
tion point Ec ~ 3.4 the spectra are qualitatively similar 
to the corresponding ones of the 2D model, whereas at 




FIG. 12. Local spectral density po{uj) of the 2D 
Chalker-Coddington network of 50 x 50 nodes in the local- 
ized regime a,t E = 2Et (upper) and at the transition point 
(lower) aX E = Ec = Q. The localized density is concentrated 
on few frequencies, indicating strong localization. The critical 
density is distributed over the entire band and appears to be 
scale invariant. 
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FIG. 13. Local spectral densities of a 3D network model of 
50^ nodes in the localized {E — QEt), critical {E — 3.4) and 
metallic {E — 1) regime. 

energy E = 1, where states are delocalized, the spectrum 
is comparatively homogeneous. 

The critical spectral densities deserve a closer look. 
Their scale invariance (or self-similarity) can be demon- 
strated by examining the scaling behavior of the quan- 
tity m{Q) = J " pq(lu) dcu, which is the mass of the 
interval [ujqjLOo + Q] with respect to the density pa{i-j). 
Disorder averaged powers of it, mq{n) ~ (m'^(r2)), scale 
over a wide frequency (Q) range with definite exponents, 
as can be seen by the linear dependence of himq(Q) on 
In Q, plotted in Fig. |lj. This power law scaling of ■m'^{fl) 
proofs the scale invariance of the density po(^)- Unlike 
in the case of homogeneous densities, the scaling expo- 
nent depends nonlinearly on q, therefore the densities are 
called multifractals [Q. They can be characterized by 
generalized dimensions Dq describing the scaling of the 
mq{Q) by mq{Q) oc Qd.+(i-^')^'i ^ for, in principle, all real 
q. d denotes the dimension of the density's support, thus 
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FIG. 14. Scaling behavior of qth powers of masses 
m{Q,) = f^" po{^) duj for box sizes Q. from 27r/4096 up 
to 2-n and powers q ranging from -2 to 3 in steps 0.5. The 
data result from the average over non overlapping intervals of 
28 systems of 200 x 200 nodes each. 



d = 1 for the spectral measure. [| Figure |15| shows the 
generalized dimensions of the critical spectral densities 
of both, the 2d and the 3d model. For the determination 
of the Dq we use 28 2d systems of 200 x 200 nodes and 
5 3d systems of 50 x 50 x 50 nodes, respectively. The 




FIG. 15. Generahzed dimensions of critical local spectral 
densities in the 2D (O) and 3D (D) network model. 

values for some characteristic exponents, namely the in- 
formation dimension Di, the correlation dimension Z32, 
D^oo, Doo, and the quantity ao == ^ ~ dDq/dq{q = 0) 
are listed in Tab. |. The coincidence of the correlation 



* Equivalently, instead of the generalized dimension Dq often 
also exponents r, = (g — 1)^9 or the Legendre transformed 



Dim. 


Di 


D2 


D-cc 


D^ 


ao 


2rf 


0.87 ±0.01 


0.75 ±0.01 


1.5 ±0.2 


0.6 ±0.1 


1.14 ±0.01 


3d 


0.67 ±0.02 


0.48 ± 0.03 


2.1 ±0.1 


0.2 ±0.2 


1.34 ±0.2 



/(a) of the function r, 
density H. 



are used to describe a multifractal 



TABLE I. Generahzed dimensions and ao of critical local 
spectral densities in the 2D and 3D network model. 



dimension D2 = 0.75 ± 0.02 and 0.48 ± 0.03 respectively 
with the decay exponents of the return probability p{t), 
X = 0.76 ± 0.01 and 0.43 ± 0.04 is not by accident, but 
because Po[ijj) and p{t) = | -00(^)1^ are related via V'o(i) 
by Fourier transformation, which implies that x equals 



B. Quasi-Energy Level Distribution in the 
Chalker-Coddington Network 

By interpreting C/ as a time evolution operator, the 
eigenphases loi of the unimodular eigenvalues e"^' of U 
must be considered as the analog to energy levels. In 
the following we will consider the distribution of these 
eigenphases or quasi-energy levels LOi{e) of the Chalker- 
Coddington network. We will determine the level spacing 
distribution function P{s) and the level number variance 
Yi2{N). P{s)ds gives the probability to find a quasi- 
energy level at distance sds of another level, whereby 
s is measured in units of the average level spacing A. 
Yj2{N) = {{n — (n))^) denotes the variance of the num- 
ber of levels in an interval which contains on average 
(n) = N levels. 

The investigated systems are closed networks with pe- 
riodic boundaries oi LxL — 50x50 nodes in the localized 
and critical regime at energies E — 10 and E = Ec = 0, 
respectively. The corresponding unitary network opera- 
tors U{E) are represented by A/ = 2 x 50 x 50 dimensional 
matrices. Diagonalizing such a matrix by standard nu- 
merical methods yields M complex eigenvalues 6*"^'*^^^ 
homogeneously distributed on the unit circle. The cor- 
responding quasi-energy levels ujn (E) are homogeneously 
distributed in the interval [0, 27r[ with an averaged spac- 
ing A — 2tt/M. It is not necessary to perform spe- 
cial unfolding-procedures in order to compensate inho- 
mogeneities in the distribution. For these quasi-level dis- 
tributions uji{E) we determine P{s) and £2(5). 

Fig. |l^ shows the level distribution function P{s) for 
networks in the localized regime at energy ii^ = 10 and at 
the transition point at _E = £'c = 0. Since for £' = 10 the 
networks states are strongly localized, their eigenvalues 
(quasi-energies) are uncorrelated, reflected by the Pois- 
son distribution P{s) — cxp(— s). Accordingly, the level 
number variance is given by T,2{N) — N (Fig. p7| ). 

At the transition point, eigenstates are extended giving 
rise to a strong level repulsion, as can be seen by the cor- 
responding level distribution function, which resembles 
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FIG. 16. The level distribution function P{s) for 
quasi-levels of networks in the localized regime, E — 10 (x), 
and at the transition point E = Ec = (+)• In the localized 
regime the distribution is clearly governed by Poisson statis- 
tics (solid curve), whereas it is close to the Wigner surmise 
for a unitary ensemble (dashed line) at the transition point. 



the Wigner surmise for a unitary ensemble (Fig. 16). 
However, Fig [l/^ shows that the level number variance is 
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FIG. 17. Level number variance E2(A'^) for quasi-levels 
of networks in the localized regime, E — 10 (x), and at the 
transition point, E — {+). 

stih linear in N, E2(^) = (0.124 ±0.006) x N. This is in 
contrast to metallic behavior, where S2(A^) oc In(A^), but 
resembles the localized phase where E2(iV) = N. The 
difference to the latter is only the reduced proportional- 
ity constant, the spectral compressibility, x — 0.124 < 1. 
This value of the spectral compressibility fits to the result 
of Chalker et al. pj], according to which x — (2 — D2)/2d, 
D2 being the correlation dimension of critical states in 
Quantum Hall systems uM- The non- vanishing com- 



pressibility is reflected by a longer tail in the level spacing 
distribution. As predicted by Altshuler et al. [^ the tail 
of P(s) decreases as exp(— ks), where k = l/(2x) ~ 4. In 
p3| , p^ 15 1 this relation could be verified for the Anderson 
model and the Chalker-Coddington model, respectively. 
Since there is no real metallic phase in the QHE, we can- 
not give examples for metallic quasi-level distributions 
from the Chalker-Coddington network. 



For the 3D network this is possible. Fig 18 shows that 



the levels spacing distribution in the localized and metal- 
lic regime at energies E = and E = 10 agree well 
with Poisson statistics and the Wigner surmise for a uni- 
tary ensemble, respectively. (Data obtained from 40000 
quasi-energy levels of 10 x 10 x 10 nodes networks.) The 




FIG. 18. Level distribution function P{s) for quasi-levels 
of 3d networks in the localized {E = 10, x ) and metallic 
{E — 0,+) regime. The data fit well to Poisson statistics 
(solid curve) and to the Wigner surmise for a unitary ensemble 
(dashed line), respectively. 

corresponding level number variances plotted in Fig. 03 
indicates Poisson statistics in the localized regime and 
show the strongly suppressed variances for the metallic 
regime. 

These two examples demonstrate that the distribution 
of network specific quasi-energy levels is equivalent to the 
distribution of real energy levels, as could be anticipated 
from the fact that network models correctly reproduce 
the dynamical properties of real physical systems. 

Another argument supporting this equivalence relies 
on the symmetries and homogeneity of the fiow of the 
quasi-energy levels lui (E) as a function of the real energy 
E. It can be found in p^ . 

VI. EIGENSTATES OF THE NETWORK 

Correlations between eigenvalues of random matrices 
come along with certain correlations in the amplitudes 
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FIG. 19. Level number variance E2(A'^) for quasi-levels of 
3d networks in the localized ( x ) and the metallic (+) regime. 



of eigenstates, which will be the subject of this section. 
As before, we investigate correlations of the Chalker- 
Coddington network, whereby we confine ourself to the 
critical region. Our findings for the network model will 
be compared to results obtained within other, conven- 
tional models, describing Quantum Hall systems with 
short range disorder by a tight-binding Hamiltonian pq ] 
and by a so-called random Landau model |g^ . 

Network eigenstates ^„ at energy E are solutions of 
the eigenvalue problem 



U{E)^„ = e''^"*„. 



(30) 



Notice that the energy E is fixed and only the quasi- 
energy Un varies. A particular eigenstate '^„ with link 
amplitudes {ipi}i=i^...N refers to a spatial wave function 
\l/(r) via the link locations r; by identifying ^(n) = ^pl. 
This rather vague definition is satisfactory as long as 
the considered length scales are large compared to mi- 
croscopic scales. 

Fig. po| shows the square-amplitude p{r) — |^(r)p of 
a typical eigenstate obtained in the Chalker-Coddington 
network at critical energy E = Er = 0. Like the lo- 



cal spectral density (see Sec. |VA| ), the eigenstate ex- 
hibits a self-similar, multifractal structure, generic for 
the localization-delocalization transition point. Similar 



FIG. 20. Squared amplitude of a Critical wave functions 
in a Chalker-Coddington network of 256x256 saddle points. 
Darker areas denote lower square amplitude. 



as for the local spectral density, the scale invariance of 
the eigenf unctions can be demonstrated by investigating 
the scaling behavior of the quantity 
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FIG. 21. Power-law scaling of the average qth-moment of 
box masses mjZ) for g-values from -2.0 to 3.0. 



ra{l) 



(frp{ro +r), 



(31) 



which is the mass of a square of linear dimension / with 
respect to the density p{r). The disorder averaged qth- 
moments of it, mq{l) = {mf{l)), obeys power law behav- 
ior over a wide range of box-sizes ?, as can be seen in 
the double logarithmic plot in Fig. |2l|. T he data were 
taken from the wave function in Fig. |20[. The general- 
ized dimensions Dq of the critical density p(r), defined by 
■mq{l) ex Zd-f(?-i)i5, (J ^ 2), are plotted in Fig. H. For 
the determination of these exponents we used networks 
with sizes varying from 80 x 80 up to 512 x 512 and of two 
different kinds of boundary conditions, full-periodic and 
semi-periodic ones. To calculate eigenfunctions we ap- 
plied the method of inverse iteration and, to treat larger 
systems, the Lanczos algorithm. Significant deviations in 
the multifractal exponents for networks of different sizes 
and boundary condition were not found. The results are 
also independent on the used numerical method. Table O 
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22. The function D{q) averaged over 12 wave func- 
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lists some characteristic exponents and gives a compari- 
son to the results obtained in the two models mentioned 
above. The exponent ao is defined as the derivative of Dq 
with respect to g at g = 0. The values for the multifrac- 
tal exponents of these models match within the errors. 
The corresponding exponents of the critical local density 



Model 



N 



TB 



RL 



D/D 



Di 



1.75 ±0.01 



2.01 ±0.02 



D2 



1.50 ±0.01 



1.62 ±0.02 



1.5 ±0.1 



2.00 ±0.02 



D-c 



3.4 ±0.2 



3.7 ±0.2 



3.7 ±0.02 



2.3 ±0.4 



Do 



1.1 ±0.1 



1.0 ±0.1 



0.95 ±0.2 



1.8 ±0.4 



ao 



2.26 ±0.01 



2.29 ±0.02 



2.3 ±0.07 



1. 



1±0.02 



TABLE II. Comparison of characteristic multifractal expo- 
nents for different models: network (N), tight-binding (TB) 
and the Random Landau (RL) model . The last row gives the 
ratio of corresponding exponents for critical eigenstates and 
critical local spectral density. 

of states (with tilde) appear to be half of those of the 



9' 



critical states. This simple proportionality, Dq = dD, 
where d = 2 is the spatial dimension, is a consequence of 
dynamical scaling introducing a single energy dependent 
length scale L^ oc w"^/'' Q. 



VII. TWO-POINT-CONDUCTIVITIES 

Finally, we come to an example of an open network. 
We will study a Chalker-Coddington network coupled to 
external reservoirs via ohmic contacts, in order to demon- 
strate how the network model can be used for the deter- 
mination of conductivities. For simplicity we assume the 
the coupling to the reservoirs is given by two ideal point 
contacts: two links, a and &, of a closed network in a 
distance r are cut and connected to external reservoirs of 
chemical potential y.ailJ'b, respectively see Fig. g^. This 
means that, firstly, the incoming channels a' and h' are 
occupied up to the chemical potentials ^ia,l^b^ secondly, 
the currents in the outgoing channels a and h are ide- 
ally absorbed by the reservoirs. The goal is to calculate 
the (dimcnsionless) conductance g of this configuration 
obeying the equation 51 = e/h gS^, where S^j. = fia — ^J■b 
and SI the net current flowing through the network from 
reservoir a to b. Here the entire network can be viewed 
as a single scattering center - with quite complex in- 
ternal structure - connecting the two incoming channels 
{a',b') to the outgoing ones (a, 6). So, according to the 
Landauer-Biittiker formula, the conductance g is in this 
case the corresponding scattering amplitude \tab\'^ from 
a' to 5, which to determine is our task. To do this, con- 
sider the following situation: let the incoming channel a' 
be feed by a unit current with amplitude 1 and set the 
amplitude in the other one, b' , to 0. The current will dis- 
tribute over the network obeying (i) the boundary condi- 
tion tpa' = 1, ^b' = at the external incoming links a' ,b' 




FIG. 23. Chalker-Coddington network coupled to external 
reservoirs by two point contacts. 



and (ii) the local scattering conditions f/;/ = '^„^ timtpm 
at internal links I (including a and b). Then, ipb, the 
amplitude in outgoing channel b, equals exactly tab, the 
coefficient we are looking for. It is straightforward to 
show that conditions (i) and (ii) are equivalent to the 
vector equation 



C/(* - t/jaea + ea - ipbeb) = *, 



(32) 



where all quantities refer to the closed network. The term 
—4'aGa belongs to the current reflected back into reservoir 
a, +ea corresponds to current fed into channel a', and 
—il^b^b is due to current out-flowing through channel b' . 
With projection operators P/ = 1 — P; = 1 — |e/ >< e;| 
Eq. (^) reduces to 

(1 ~ UPaPb)^ = C/e„ 
and we obtain finally, 

tab = V'h = (efc,*) 

^{eb,{l-UPaPby^Uea). 

By this, the problem of finding the two-point conduc- 
tance is reduced to a standard problem of linear algebra, 
which we solved by numerical LU-decomposition. In the 
following we will present data obtained by this method 
for critical Chalker-Coddington networks. 

At criticality the conductances are expected to show a 
power law behavior with respect to the distance r <^ L 
between the contacts, while their distribution will be very 
broad for any r. Furthermore, as was pointed out by 
Zirnbauer H, there should not be a dependence of the 
conductance on the system size. As we will see, such a 
picture is consistent with our data. We have determined 
conductances for systems of size L — 40, 60, 100 and dis- 
tances between links varying from ?' = 1 to r = L/2. For 
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every distance we determined up to 1300 conductances 
(depending on the system size). 

The broadness of the distribution suggests to deter- 
mine the moments [g''], the geometric mean gt — exp[ln 5] 
and the log- variance [(Ing/gt)^] . In Fig. e3 we show 
[In 17] plotted against Inr for three different system sizes 
L — 40, 60, 100. There is no visible L dependence, as ex- 
pected, and the typical value scales in a power law fashion 
with respect to r. 
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FIG. 24. The logarithm of the typical conductance plotted 
against In r for three different system sizes. 



We find similar behavior for the moments [g'], whose 
behavior can be described in the following way: 



[g'^]L oc r 



-X(q) 



(33) 



where X{q) is called the multifractal spectrum of the 
conductance. The numerical data for X{q) for system 
sizes 40, 60 and 100 is depicted in Fig. pM Within the 
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FIG. 25. The function X{q) for three different system sizes. 



error bars, the non-linear curves for different system sizes 
match, indicating universal multifractality. For further 



information on the distribution P{g) of the point-contact 
conductances see ref. BSl. 



VIII. SUMMARY 

We extended the Chalker-Coddington model and a 3D 
generalization of it by a simple, discrete time evolution 
and studied dynamical and spectral aspects of the result- 
ing dynamical network models. Despite of the simplicity 
of the dynamics, the models correctly show the univer- 
sal aspects of normal and critical diffusion as well as of 
Anderson localization in temporal evolution of network 
wave functions. In case of closed networks, the analogue 
of eigenf unctions and eigenenergies exist in form of eigen- 
vectors and quasi-energies (=eigenphases) of an unitary 
network operator. We investigated spatial and spectral 
correlation of eigenstates and quasi-energy spectra and 
found qualitative and quantitative agreement with theo- 
retical predictions and other numerical work within con- 
ventional models. Furthermore, we used an open network 
for the calculation of two-point conductances. 

Definitions and methods with regard to the network 
dynamics are not specific to the two models investigated 
in particular. Therefore, they should be as well applica- 
ble to any other model formulated as a network in terms 
of local scatterers and connecting link-channels. 

Disregarding physical and mathematical subtleties, the 
analogies of a physical system described by a Hamiltonian 
H with it's corresponding network model represented by 
an unitary operator U can be summarized as in Tab. 
III. Notice that there are still significant differences: 



g-iHT 


u 


*(f) = e-"^'*(0) 


^{kr) = [/'=*(0) 


Htpn = £„1p„ 


(/^^^e''^"*^ 


. . . £71 — l,£n,£n + l, . . . 


. . . OJn-l, l-^n, (^n + 1, ■ ■ ■ 



TABLE III. Relation between Hamiltonian and network 
operator (schematic). 

first of all, the dimensions of H and U do not match. 
The dimension of U equals the number of links, which 
generally is much smaller than the dimension of if, since 
each link-channel itself contains a large number of (quasi- 
) asymptotic states. In some sense, changing from H to 
U corresponds to reducing the original system by non- 
interesting degrees of freedom. In case of the models 
investigated here, these degrees of freedom are the link 
states: the network wave function \1/ no longer contains 
any information about their structure (e.g. wave num- 
ber). Obviously, this dimensional reduction constitutes 
the main advantage of the network model, especially in 
numerical applications. The price to be paid for that is a 
loss of accuracy. Processes and correlations correspond- 
ing to length and energy scales of the out-projected in- 
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ternal degrees of freedom can no longer be reproduced. 
In particular this means that the statistical equivalence 
of energy and quasi-energy spectrum holds only in the 
vicinity of a certain energy Eq to which U = U{Eo) refers. 
Similarly, the network time evolution can only be related 
to the time evolution of real wave packets composed of 
eigenstates with energies close to Eq. 
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